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Abstract 


Two  major  areas  of  effort  are  encompassed: 

I,  New  Techniques  for  the  Synthesis  of  Metals  and  Alloys 

The  High  Rate  Physical  Vapor  Deposition  (HRPVD)  process  is  to  be 
used  for  the  following: 

1.  Preparation  and  characterization  of  Ni  and  Ni-20Cr  alloy  sheet 

2.  Synthesis  of  compounds  Y2®3*  ^iC,  S*3N4  by  reactive  evaporation 
and  their  characterization 

3.  Dispersion  strengthened  alloys,  Ni-20Cr-Y20^,  Ni-20Cr-TiC  and 
Ti-Y203 

This  report  describes  the  HRPVD  apparatus.  Models  for  calculation  of  tempera¬ 
ture  and  thickness  distribution  of  the  deposit  are  given. 

II.  The  Properties  of  Rare  Earth  Metals  and  Alloys 

Small  amounts  of  certain  rare-earth  elements  in  Nickel-base  alloys 
are  known  to  markedly  reduce  the  rate  of  oxidation  at  high  temperatures  and 
confer  mechanical  stability  on  the  scales  found.  This  program  is  directed 
at  obtaining  an  understanding  of  the  mechanisms  and  to  characterize  the 
oxidation  behavior  of  Ni-20Cr  containing  Gd,  La  and  Y  and  alloys  based  on 
composition  Ni^Al. 
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1.  Introduction 

This  report  describes  research  activities  on  ARP  A  Grant  No.  A  01643. 

The  scope  of  the  woik  is  divided  into  two  major  areas  of  effort  and  further  sub¬ 
divided  into  four  tasks  as  shown  below. 

1.  New  Techniques  for  the  Synthesis  of  Metals  and  Alloys  -  Tasks  I,  II, 
and  III.  (Professor  R.  F.  Bunahah  -  Principal  Investigator) 

2.  The  Properties  of  Rare  Earth  Metals  and  Alloys  -  Task  IV.  (Professor 
D.  L.  Douglass  -  Principal  Investigator) 

In  the  following  each  of  the  two  areas  and  tasks  will  be  described 
separately  with  the  progress  to  date. 

II.  New  Techniques  for  the  Synthesis  of  Metals  and  Alloys  -  (Tasks  I,  II 
and  III) 

A.  Background 

High  Rate  Physical  Vapor  Deposition  (HRPVD)  Techniques^1 ^ 
are  to  be  used  to  prepare  metallic  alloys,  ceramics,  and  metal-ceramic 
mixtures  (dispersion  strengthened  alloys) .  The  method  consist  s  of  evapora¬ 
tion  of  metals,  alloys  and  ceramics  contained  in  water  cooled  crucibles  using 
high  power  electron  beams.  The  process  is  carried  out  in  a  high  vacuum 
environment.  The  use  of  high  power  electron  beams  makes  it  possible  to 
produce  very  high  evaporation  rates.  The  vapors  are  collected  on  heated 
metallic  substrates  to  produce  full  density  deposits  at  high  deposition 
rates. 

There  are  three  (3)  tasks  in  this  section: 

Task  I:  The  preparation  and  characterization  of  Nickel  and  Ni-20Cr 

alloy  sheet  by  the  high  rate  physical  vapor  deposition  process. 

Task  II:  Synthesis  and  characterization  of  compounds  by  Reactive  Evapora¬ 

tion.  The  compounds  to  be  prepared  are  YjO^,  TiC  and  Si^N^. 

Task  III:  Dispersion  strengthened  alloys  produced  by  HRPVD  Process,  and 

their  characterization.  The  specific  alloys  to  be  studied  are: 

1.  Ni-20Cr-Y203  2.  Ni-20Cr-TiC  3.  Ti-Y^ 
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Single  source  and  two  source  evaporation  methods  will  be  used 
to  produce  these  alloys. 

The  HRPVD  process  has  several  attractive  features: 

1.  Simple,  full  density  shapes  (sheet,  foil,  tubing)  can  be  produced  at 
high  deposition  rates ,  5s- 0.001"  per  minute  thickness  increment  thus  making 

f 

it  an  economically  viable  process. 

2.  Metals  and  alloys  of  high  purity  can  be  produced. 

3.  Very  tine  grain  sizes  (1/U grain  dia.  or  smaller)  can  be  produced  by 
controlling  substrate  temperature.  Grain  size  refinement  is  produced  by 
lowering  the  condensation  temperature. 

4.  An  alloy  deposit  may  be  produced  from  a  single  rod  fed  source.  This 
occurs  because  the  molten  pool  at  the  top  of  the  rod  is  only  1/4"  deep. 

The  vapor  composition  is  the  same  as  that  of  the  solid  rod  being  fed  into 
the  molten  pool.  At  equilibrium,  the  composition  of  the  molten  pool  differs 
from  that  of  the  vapor  or  the  solid  feed.  It  is  richer  in  those  components 
having  a  low  vapor  pressure.  The  composition  of  the  vapor  is  the  product  of 
the  vapor  pressure  times  the  mol  fraction  of  the  component.  For  example, 

a  T1-6A1-4V  alloy  deposit  where  the  differences  in  vapor  pressure  of  A1  and 
V  are  a  factor  of  5,000  at  1600°C  can  be  produced  by  evaporation  from  a 
single  source.  The  feed-rod  is  T1-6A1-4V  and  the  molten  pool  is  much  richer 
in  V  than  in  Al. 

5.  Two  or  more  sources  can  be  used  to  simultaneously  deposit  on  the  same 
substrate  thus  conferring  the  ability  to  produce  complex  alloys.  For  example, 
an  alloy  with  a  2  or  3  component  solid  solution  matrix  may  be  evaporated 
from  one  source  and  another  metal  or  ceramic  for  the  dispersed  phase  from 
another  source.  The  dispersion  size  and  spacing  should  be  very  fine  since 
the  deposition  is  occurring  from  the  vapor  phase. 

The  unique  feature  of  this  process  is  that  all  of  the  above  benefits  can  be 
obtained  simultaneously. 


It  should  be  noted  that  the  condensation  temperature  is  a  very  impor- 

(8) 

tant  process  variable.  Bunshah  and  Juntz'  '  found  that  for  titanium^ as  the 
deposition  temperature  is  lowered  the  grain  size  of  the  fully  dense  deposit 
becomes  finer.  At  very  low  temperatures  (fV25%  of  the  melting  point)  the 
deposit  has  less  than  full  density.  Since  a  fine  grain  sized  microstructure 
represents  an  optimum  condition  of  strength  and  toughness  in  a  material,  the 
Importance  of  control  of  the  deposition  temperature  becomes  obvious. 

B.  Scope  of  Work 

The  main  tasks  on  this  contract  are  the  preparation  and  testing  of 
the  various  alloys,  ceramics  and  dispersion  strengthened  alloys  as  outlined 
in  Section  I  above.  Very  essential  to  the  preparation  of  suitable  test 
specimens  are  two  other  factors: 

1.  Design  of  the  apparatus  for  high  rate  physical  vapor  deposition. 

2.  Theoretical  calculation  of  the  thickness  distribution  and  temperature 
distribution  of  the  deposited  material  which  in  this  case  is  in  the 
form  of  a  sheet. 

Both  of  these  tasks  are  essential  preliminaries  to  the  main  scope  of 
work  and  are  detailed  below. 

C.  High  Rate  Physical  Vapor  Deposition  Apparatus  for  Alloys, 
Dispersion  Strengthened  Alloys  and  Ceramics 

Such  an  apparatus  is  quite  complex  and  has  many  state-of-the-art 

features.  The  design  requirements  for  such  an  apparatus  are  as  follows: 

1.  High  rate  evaporation  of  metals,  alloys  or  ceramics  from  one  source  or 
two  sources  simultaneously.  Sources  should  be  rod-fed  to  permit  long  time 
uniform  evaporations. 

2.  Electron  beam  heating  incorporated  into  the  sources  to  get  high 
evaporation  rates.  In  addition,  beam  scanning  of  the  source  should  be 
available  in  one  of  the  sources  which  is  used  for  the  direct  evaporation 


of  ceramics  or  dielectrics 
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3.  Substrate  holder  capable  of: 

a)  Heating  the  substrate  to  an  elevated  temperature  and  controlling  It 
at  the  desired  temperature.  This  is  necessary  to  produce  a  fully 
dense  deposit.  The  apparatus  will  have  capability  of:- 

i.  heating  a  10"  x  10"  substrate  up  to  1000°C  by  radiant  heating 
from  the  bacx. 

ii.  direct  resistance  heating  of  one  or  more  1/2"  x  3"  strips  to 
1500°C. 

b)  Rotation  of  the  10"  x  10"  substrate  and  the  ability  to  tilt  the 
axis  of  rotation  up  to  45°  from  the  vertical.  This  is  necessary 

to  produce  a  uniform  thickness  of  deposit  over  the  entire  substrate. 

4.  A  high  vacuum  environment  capable  of  maintaining  a  working  pressure  of 
5.10  ^  tcrr  or  lower  and  a  base  pressure  of  2.10  ^  torr  or  lower.  The 
vacuum  environment  should  be  monitored  using  residual  gas  analyzers. 

The  best  vacuum  techniques  should  be  used  to  keep  contamination  of  the 
deposit  from  back-streaming  and  other  sources  to  a  minimum. 

5.  A  liquid  pool  level  monitor  to  keep  the  height  of  the  liquid  pool  at 
the  desired  level.  This  is  necessary  to  keep  the  alloy  composition  in 
the  vapor  phase  at  a  constant  value. 

6.  Capability  to  bleed  gases  or  gas  mixtures  at  the  desired  partial  pressures 
in  the  vicinity  of  the  substrates.  This  is  necessary  to  carry  out 
reactive  evaporation  processes  effectively. 

7.  The  ability  to  feed  raw  material  and  melt  it  into  evaporant  stock.  This 
feature  is  very  desirable  as  the  raw  material  for  evaporant  stock  should 
have  very  low  gas  content  to  avoid  spattering  of  liquid  droplets  during 
evaporation.  Not  all  evaporant  materials  can  be  purchased  in  the  desired 
form.  Rod  fed  electron  beam  sources  also  require  the  feed  rod  to  have 
outside  diameter  to  fit  the  inside  diameter  of  the  opening  in  the  source 
after  allowing  for  thermal  expansion.  Evaporant  materials  may  not  be 
available  in  the  desired  rod  sizes. 

8.  The  ability  to  measure  the  temperature  of  the  molten  pool  preferably 
with  an  optical  or  total  radiation  pyrometer.  This  requires  clear 
visual  access  to  the  molten  pool  and  the  ability  to  keep  the  window 
free  of  deposited  vapors.  The  latter  devices  are  commercially 
available. 
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D.  Temperature  Distribution  and  Thickness  Distribution  of  the  Deposit 

Consider  a  flat  substrate  located  at  some  distance  from  an  evaporation 
source.  For  a  point  source  or  a  small  flat  plate  source  the  vapor  density 
will  decrease  proportional  to  cos  0,  where  8  Is  the  angle  between  the 
vertical  axis  of  the  source  and  the  location  on  the  substrate  (Cosine 
distribution  law) .  Thus  the  deposit  will  be  tapered  from  the  center  to 
the  edge.  This  is  not  too  desirable  as  it  will  affect  the  specimen  yield 
(i.e. ,  the  number  of  specimens  that  can  be  machined  from  a  deposit  without 
grinding  it  to  uniform  thickness)  as  well  as  the  temperature  distribution 
from  the  center  to  the  edge  of  the  specimen.  The  percent  thickness  variation 
from  center  to  edge  on  a  deposit  varies  with  the  distance  between  source  and 
substrate  and  usually  becomes  less  than  10%  at  source-substrate  distance  of 
10  times  the  source  diameter.  In  electron  beam  heated  sources  the  thickness 
variation  is  somewhat  greater  than  that  produced  by  other  sources  and  also 
depends  on  the  particular  material  being  evaporated.  The  total  thickness 
per  unit  time  also  varies  inversely  as  the  square  of  the  source-substrate 
distance.  Therefore  a  calculated  estimate  of  the  deposit  thickness  and 
thickness  variation  as  a  function  of  source-substrate  distance  is  necessary 
for  proper  planning  of  the  experiment.  One  method  of  producing  uniform 
thickness  deposits  is  to  collect  them  on  rotating  substrates  with  the  rota¬ 
tion  axis  inclined  about  30°  to  the  horizontal  plane.  For  electron  beam 
sources,  the  exact  angle  would  have  to  be  determined  by  experimentation. 

The  deposition  temperature  is  very  important.  At  very  low  deposition 
temperatures  the  deposit  will  not  be  fully  dense.  ut  higher  deposition 
temperatures,  the  deposit  will  be  fully  dense.  The  grain  size  of  the 
deposit  will  coarsen  as  the  deposition  temperature  gets  higher.  Since  the 
strength  and  toughness  of  a  deposit  decrease  with  increasing  grain  size, 
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the  selection  of  the  deposition  temperature  becomes  a  very  important  process 
variable.  It  should  be  high  enough  to  produce  a  fully  dense  deposit  and 
yet  not  too  high  to  permit  grain  size  coarsening. 

An  estimation  of  the  temperature  of  a  deposit  produced  by  high  rate 
physical  vapor  deposition  processes  is  somewhat  complicated.  A  new  method 
has  been  developed  for  making  such  a  calculation.  By  performing  an  energy 
balance  between  energy  gain  and  energy  loss  of  the  deposit.  The  deposit 
receiver  energy  from  three  sources  - 

1.  Radiation  from  the  hot  source. 

2.  Latent  heat  of  condensation  of  the  vapor  and  cooling  of  the 
deposit  to  the  deposition  temperature. 

3.  Heat  of  reaction  when  a  deposit  species  reacts  with  the 

gas  phase  to  form  a  compound.  This  is  applicable  in  the  case 
of  reactive  evaporation  only. 

In  a  vacuum  environment  where  the  substrate  is  supported  so  as 
to  minimize  conductive  heat  transfer  to  the  holder,  the  principal  mode  of 
heat  loss  is  radiative.  Radiative  heat  loss  occurs  from  the  deposit  face  as 
well  as  the  back  side  of  the  substrate.  In  addition,  there  is  a  small  energy 
loss  necessary  to  heat  the  substrate  to  the  experimental  temperature;  however 
this  factor  is  small  for  these  substrates  and  is  zero  once  equilibrium  has  been 
established.  The  time  to  attain  equilibrium  is  usually  less  than  a  minute  and 
in  most  cases  a  very  small  fraction  of  the  total  deposition  time. 

The  experimental  variables  are: 

1.  Nature  of  evaporating  species  (metal,  alloy,  compound)  and  the  vapor 
pressure  of  the  component  species. 

2.  Temperature  of  the  source  and  evaporation  rate. 

3.  Source  diameter. 

4.  Source-substrate  distance. 

5.  Substrate  material. 

6.  Substrate  material  thickness. 
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7.  Thermodynamic  properties  of  condensate  (latent  heat  of  condensation, 
specific  heat,  vapor  pressure). 

8.  Heat  of  reaction  in  the  case  of  reactive  evaporation. 

9.  The  total  hemispherical  emittance  of  the  evaporant  as  a  function  of 
evaporation  temperature. 

10.  The  total  hemispherical  emittance  of  the  deposit  as  a  function  of 
condensation  temperature. 

11.  The  total  hemispherical  emittance  of  the  substrate  as  a  function  of 
condensation  temperature. 

12.  The  absorptance  of  the  deposit  as  a  function  of  condensation  temperature. 

With  a  large  number  of  variables,  it  is  essential  to  make  a  computational 
estimate  of  the  temperature  distribution  of  the  deposit  to  minimize  the 
number  of  experiments  performed.  A  computer  program  ha&  been  written  to 
make  such  an  estimate. 

These  factors  are  considered  in  detail  in  Supplements  1  and  2  attached 
to  this  report.  These  Supplements  are  self-contained  documents  dealing 
with  the  following  topics: 

1.  Model  for  calculating  the  deposit  temperature  in  high  rate  physical  vapor 
deposition  processes. 

2.  Temperature  and  Thickness  Distribution  on  the  Substrate  During  High  Rate 
Physical  Vapor  Deposition  of  Materials. 

They  are  to  be  submitted  for  publication. 

E.  Progress  to  date  (July  1,  1970  -  January  31,  1971) 

1.  High  Rate  Physical  Vapor  Deposition  Process  Apparatus. 

The  design  of  this  apparatus  was  completed,  specifications  written, 
quotations  obtained,  and  order  placed.  The  apparatus  will  be  delivered  in 
February  1971. 

2.  Theoretical  Calculation  of  the  Temperature  and  Thickness  Distri¬ 
bution  of  the  Deposit. 

A  new  method  has  been  developed  to  calculate  the  thickness  and  tempera¬ 
ture  distribution  of  the  deposit.  The  importance  of  this  has  been  explained 
above.  Supplements  1  and  2  attached  to  this  report  describe  this  work  in 
detail.  Temperature  and  thickness  variation  for  yttrium  and  yttria  O^O-j) 
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deposits  have  been  used  to  illustrate  the  calculations.  *s  Proc^uce^ 

by  reactive  evaporation  between  metal  vapors  and  oxygen  introduced  into  the 
gas  phase.  The  range  of  variables  to  be  explored  experimentally  is  defined 
from  the  calculations  for  the  deposition  of  Y  and  They  are: 

1.  Temperature  of  the  source  -  2300  to  2600°K. 

2.  Source  to  substrate  distance  -  25  to  35  cm. 


3.  Diameter  of  source  -  1.0"  to  1.5". 


4.  Substrate  material  -  buffed  stainless  steel  and  copper. 

5.  Substrate  thickness  -  .001  to  .005". 

6.  Substrate  size  (diameter)  -  8  to  20  cm. 

3.  Deposition  of  Ni-20Cr  alloy. 

Using  a  single-source  apparatus  available  at  UCLA,  fully  dense  ductile 
deposits  of  Ni-20Cr  alloy  have  been  made  with  good  success. 

F.  Future  Work 

In  the  next  half  year  period,  the  following  work  is  scoped: 

1.  Installation  and  operation  of  the  2  source  vacuum  evaporation 


apparatus . 

2.  Calculations  of  temperature  and  thickness  distribution  of  the 

deposit  for  TiC,  Si^N^,  Ni  and  Ni-2uCr  will  be  carried  out. 

3.  Development  of  a  model  relating  the  composition  of  the  vapor 

phase  in  an  alloy  deposit  to  the  composition  of  the  molten  pool.  This  model 
will  be  developed  for  rod-fed  electron-beam  heated  sources.  The  variables 
to  be  considered  are  temperature  of  the  molten  pool,  volume  of  the  molten  pool, 
and  time  to  attain  equilibrium  composition  of  the  pool. 

4.  Feasibility  studies  on  the  production  of  ^£^3*  an<*  ^^3^4 
reactive  evaporation. 

5.  Deposition  of  Ni  and  Ni-20Cr  alloy  sheets  and  study  of  their 


structure  and  properties. 
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G.  Personnel 

The  following  personnel  have  been  working  on  the  project: 

Principal  Investigator  -  Professor  R.  F.  Bunshah  (since  July  1,  1970) 
Graduate  Students  -  Raymond  Chow  and  R.  Nimmagadda  (since  July  1,  1970) 

Post  Doctoral  Fellow  -  Dr.  A.  C.  Raghuram  (as  of  February  1,  1971) 
Technician  -  Mr.  Harold  Willey  (since  November  1,  1970) 

III.  The  Properties  of  Rare  Earth  Metals  and  Alloys  -  (Task  IV) 

A.  Introduction 

The  presence  of  small  amounts  of  certain  rare  earth  elements  in 
nickel-base  alloys  is  kiown  to  markedly  reduce  the  rate  of  oxidation  at 
high  temperatures  and  to  exert  a  highly  beneficial  effect  on  the  mechanical 
stability  of  the  scales  formed.  The  effect  is  disproportionate  to  the 
amount  added.  Generally,  a  fraction  of  a  percent,  e.g.,  approximately 
0.1  to  0.3%,  is  sufficient.  In  particular,  the  resistance  to  exfoliation 
of  the  scales  becomes  outstanding. 

The  mechanism  by  which  these  effects  occur  is  not  known  and  has  been 
subject  to  rather  w*ld  speculation.  The  program  in  progress  is  directed 
at  obtaining  an  understanding  of  the  mechanism  and  to  characterize  the 
oxidation  behavior  of  certain  alloys  containing  rare  earth  additions. 

B.  Experimental  Work 

1.  Alloy  Preparation 

Melting  and  fabrication  of  alloys  was  subcontracted  to  Stanford 
Research  Institute.  Charges  of  about  400  gms.  total  were  triple  arc-melted 
in  argon  with  an  inert  electrode  on  a  water-cooled  copper  hearth.  The 
buttons  were  flipped  over  after  each  melting  and  remelted  a  total  of  three 
times.  It  is  noteworthy  that  the  alloys  based  on  the  composition  Ni^Al 
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(approximately  13  w/o  Al)  all  cracked  upon  cooling  of  the  arc-melted  buttons. 
These  alloys  were  found  to  be  extremely  brittle  and  could  not  be  successfully 
melted,  much  less  subsequently  fabricated,  in  spite  of  numerous  attempts. 

The  other  base  material  was  Nl-20Cr  to  which  was  added  Gd,  La,  and  Y, 
Initially  In  the  amount  of  1  w/o  and  subsequently  In  the  amount  of  0.5  w/o. 
These  compositions  produced  sound  buttons  which  were  amenable  to  subsequent 
fabrication. 

The  buttons  of  the  Ni-20Cr  base  alloys  were  hot  rolled  to  0.030"  thick 
sheet  at  temperatures  In  the  range  of  1200  to  1300°C.  The  surface  oxide 
was  removed  by  sandblasting,  and  the  sheet  was  cold-rolled  to  0.020"  In 
thickness. 

Oxidation  samples  were  sheared  from  the  cold-rolled  sheet,  cleaned, 
and  annealed  for  20  minutes  at  900®C.  They  were  weighed  and  then  subjected 
to  oxidation  tests. 

2.  Oxidation  Tests 

Thermal  gravimetric  analysis  was  performed  in  a  Harrop  Unit  that 
required  about  40  minutes  to  reach  a  temperature  of  1200°C.  The  weight 
gain  Increased  rapidly  with  time  during  heat  up  commencing  at  a  temperature 
of  about  500°C.  Although  the  system  could  be  evacuated  readily,  it  was  not 
feasible  to  heat  in  vacuum  and  then  introduce  the  oxidizing  environment. 

This  problem  results  from  the  difficulty  of  establishing  the  null  point  in 
the  system  which  employs  a  linear  differential  transformer  to  record 
movement  (and  hence  weight  change)  of  the  balance  arm. 

The  rapid  Increase  of  weight  at  low  temperatures  is  most  likely  due 
to  the  Initial  formation  of  N10,  the  formation  of  this  oxide  at  low  tempera¬ 
tures  being  much  more  rapid  than  the  formation  of  other  oxides,  e.g., 

Cr2^3*  NiC^O^,  or  a  rare  earth  oxide.  This  behavior  does  not  occur  at  1200°C. 
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In  order  to  circumvent  the  above-mentioned  problem,  the  sample  was 
Introduced  quickly  into  another  furnace  at  the  same  temperature  and  held  for 
one  hour.  This  permitted  the  oxides  to  form  which  are  characteristically 
formed  at  the  temperature  of  Interest.  The  sample  was  removed,  rewelghed, 
and  placed  in  the  IGA  unit.  It  was  Ihen  possible  to  calibrate  the  unit 
and  then  heat  up  with  virtually  no  change  in  weight  taking  place  during  the 
heating  period.  The  time  at  which  the  desired  temperature  was  reached 
was  then  taken  as  60  minutes,  and  the  test  was  carried  out  for  the  desired 
time  period. 

3.  Results 

Kinetics  curves  have  been  obtained  for  Ni-20Cr  and  for  Nl-20Cr-lGd 
using  the  technique  described.  Runs  were  made  at  1200°C  In  air.  The  alloy 
containing  Gd  showed  a  very  rapid  Initial  Increase  In  weight  followed  by  a 
continuously  decreasing  rate  (less  than  parabolic) .  The  presence  of  Gd 
resulted  In  a  decreased  rate  compared  to  plain  Nl-20Cr. 

C.  Personnel 

The  following  personnel  have  been  working  on  Part  II  of  this  program: 
Principal  Investigator  -  Professor  D.  L.  Douglass 
Graduate  Student  -  J.  R.  Kuenzly 


References 


1.  Bunshah,  R.  F. ,  "Superpurification  of  Metals  by  Vacuum  Distillation: 

A  Theoretical  Study,"  Trans.  Vac.  Met.  Conf.,  1963,  AVS,  121. 

2.  Bunshah,  R.  F. ,  and  Juntz,  R.  S.,  "Purification  of  Beryllium  by 
Crucible  Free  Melting  and  Distillation  Process"  in  Beryllium  Technology, 
Gordon  and  Breach,  1964,  1. 

3.  Bunshah,  R.  F. ,  "Impurity  Removal  by  Distillation  of  Beryllium  from 

the  Solid  State,"  Proceedings,  Int'l.  Conf.  on  Beryllium,  Grenoble, 
France,  1965,  Presses  Universltarles  de  France,  108  Blvd.  St.  Germain, 
Paris,  6,  63. 

4.  Bunshah,  R.  F. ,  and  Juntz,  R.  S.,  "The  Purification  of  Beryllium  by 
Vacuum  Melting  followed  by  Vacuum  Distillation  in  an  Electron  Beam 
Furnace  with  Simultaneous  deposition  of  Sheet,"  Trans.  Vac.  Met. 

Conf.,  1966,  AVS,  209. 

5.  Bunshah,  R.  F.,  "The  Effect  of  Purification  on  Some  Mechanical 
Properties  of  Beryllium,"  Metals  Engineering  Quarterly,  Nov.  1964,  8. 

6.  Bunshah,  R.  F. ,  and  Juntz,  R.  S. ,  "Electron  Beam  Distillation  Furnace 
for  Reactive  Metals:  Design  Considerations  and  Operating  Experience," 
Trans.  Vac.  Met.  Conf.,  1965,  AVS,  200. 

7.  Bunshah,  R.  F.  and  Juntz,  R.  S.,  "Design  Considerations  for  the 
Production  of  Massive  Deposits  of  Alloys  by  Evaporation  from  Multiple 
Electron  Beam  Heated  Sources,"  Trans.  Vac.  Met.  Conf.,  1967,  AVS,  799. 

8.  Bunshah,  R.  F.,  and  Juntz,  R.  S.,  "Influence  of  Condensation  Temperature 
on  Microstructure  and  Mechanical  Properties  of  Titanium  Sheet,"  to  be 
published. 


1.1 


SUPPLEMENT-1 


Model  for  Calculating  the  Deposit  Temperature  in 
High  Rate  Physical  Vapor  Deposition  Process 


Objective 

The  high  rate  physical  vapor  deposition  process  is  used  to  produce 
metals  or  alloys  of  high  purity,  full  density,  small  grain-size  and 
controlled  thickness.  The  substrate  temperature  during  deposition  is 
the  most  important  experimental  parameter  for  controlling  the  grain-size 
and  density  of  the  deposit.  As  grain-size  decreases,  the  strength  and 
toughness  of  the  material  increase.  Therefore,  it  is  important  to 
develop  a  model  for  calculating  the  temperature  distribution  of  the  de¬ 
posit.  In  this  paper,  we  will  postulate  a  model  for  calculating  the 
temperature  of  the  deposit  and  then  proceed  to  verify  the  model  by 
calculating  the  deposit  temperature  for  experimental  conditions  given 
in  the  literature  for  the  evaporation  of  titanium^  and  comparing  the 
calculated  temperature  to  the  measured  temperature.^  The  temperatures 
of  the  deposit  and  the  substrate  on  which  it  is  deposited  are  the  same 
since  they  are  in  intimate  thermal  contact.  This  has  been  experimentally 
verified  by  Bunshah  and  Juntz.^ 

Equation  of  Energy  Balance 

The  experimental  set-up  for  high  rate  physical  vapor  deposition 
processes  (HRPVD)  is  shown  in  Pig.  1.  The  metal  is  contained  in  a  water- 
cooled  copper  crucible  and  heated  by  an  electron  beam  source.  The 
metal  vapors  are  condensed  on  a  metallic  substrate  located  above  the 


source. 
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Figure  1.  Schematic  of  High  Rate  Physical  Vapor  Process 
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The  model  for  calculation  of  the  deposit  temperature  Is  based  on 
the  energy  balance  between  the  energy  Input  and  output  of  the  substrate. 

The  energy  Input  to  the  substrate  consists  of  three  terms: 

Q  -  the  radiant  energy  from  the  source  at  Temperature  T. 

K  1 

Q^,  -  energy  supplied  by  the  condensation  of  metal  vapors 

-  the  heat  of  reaction  for  reactive  evaporation  l.e.  ,  If  a 
reaction  occurs  between  the  condensing  vapors  and  a  gas 
phase  on  the  substrate  or  between  two  metal  vapors  condens¬ 
ing  on  the  substrate.  For  simple  condensation  of  metal 

or  alloy  vapors  this  term  Is  zero. 

The  energy  output  from  the  substrate  consists  of  the  following  terms: 

Rj  -  energy  radiated  from  the  deposit  side  at  temperature  T^ 
of  the  substrate 

R  -  energy  radiated  from  the  back  side  of  the  substrate  at 
temperature  l.e.,  the  side  away  from  the  source 

Hg  -  energy  necessary  to  heat  the  substrate  material  to  the 
substrate  temperature  T^ 

-  energy  necessary  to  heat  the  material  already  deposited  if  at 
a  lower  temperature,  to  the  substrate  temperature  T^ 

We  can  now  write  the  equation  of  energy  balance  as: 

«R  +  «C  +  <’te-Rd+R.  +  H.  +  Hd  (1) 

Initially  during  the  process,  the  energy  input  will  exceed  the  heat  output 
and  the  temperature  of  the  substrate-deposit  combination  will  rise,  thus 
causing  an  increased  energy  loss  from  the  substrate.  Since  the  energy  input 
to  the  substrate  deposit  combination  is  constant  and  the  output  is  propor¬ 
tional  to  T2^(T2  being  tbe  substrate  temperature),  a  steady  state  will  be 
reached.  Under  steady  state  conditions,  and  in  equation  1  are  zero 
and  the  equation  of  energy  balance  becomes 

«*  +  Rc  +  ■  Rd  +  Rs  U) 


We  are  also  assuming  here  that  the  rate  of  heat  loss  by  conduction 
of  the  substrate  and  deposit  to  adjoining  areas  is  so  small  in  comparison 
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Co  the  other  energy  loss  terms  that  It  can  be  neglected  as  shown 
In  Appendix  1. 

Let  us  now  examine  each  term  In  equation  2  In  detail. 

-  energy  absorbed  by  the  substrate  at  temperature  T2  from  the  source 
at  temperature  T 


El(X‘Tl)  W<X'V  A1  FA  -  dA  t  dX 

.  12 


(3) 


where  e^(X,T^) 


C8(X,T2) 


■  spectral  emlttance  of  source  at  temperature 

■  spectral  absorptance  at  substrate  temperature  T2 

■  Plancklan  radiant  emlttance  at  temperature  T1 


'2 

h 

k 

X 


^5  exp  (c2/  -  1 


n  2TrhC^  ■  3.7413  +  0.003  x  10  ^  watt  cm^ 

■  ^  “  1.43884  +  0.0008  cm  degree 

■  Planck's  constant 
-  Wave  number 

■  Wave  length 


A^  ■  Area  of  the  source 

dA2  -  Infinitesimal  area  at  the  center  of  deposit 


Ax  -  dA2 


ww 

t  *  duration  of  experiment 


Configuration  factor  for  radiation  traveling  from 
A,  to  dA„ 


If  we  apply  the  reciprocity  rule  for  the  configuration  factor 


(2) 


A1  FA1  -  dA2 


dA2  FdA2  -  A: 


(4) 


where  F 


dAj,  -  A^  is  the  configuration  factor  for  radiation  travelling  from 


dA2  to  A^. 
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•00 

ei(X’Tl)Qt8(X,T2)W(X*Tl)  dX 

0 


(5) 


Rd  -  energy  radiated  by  the  infinitesimal  area  dA2  of  the  deposited 


material 


Rj 


t 


*00 

-  0 


ed(X,T2)  W(X,T2)  dX 


(6) 


where  ed(X,T2) 

w(X,t2) 


spectral  emlttance  of  the  deposit  at  substrate  tempera¬ 
ture  t2 

Planckian  radiant  emlttance  at  temperature  T2> 


R  -  energy  radiated  by  the  infinitesimal  area  of  the  substrate  dA_ 


fCO 


R  *  dA0  t 
8  l 


es(X,T2)  W(X,T2)  dX 


(7) 


where 


£s^X’^2^  “  spectral  emlttance  of  the  substrate  at  substrate 

temperature  T2 


Qc  -  energy  supplied  by  the  latent  heat  of  condensation  of  metal  vapors 


Qc  *  Hs  +  h  +  Ht  (8) 

where  Hg  -  energy  released  on  cooling  of  metal  vapors  from  source 
temperature  ^  to  the  condensation  temperature  Tc  plus 
energy  released  on  cooling  the  deposited  solid  from  T 
to  substrate  temperature  T2 

-  latent  heat  of  condensation  at  temperature  Tc 
Ht  -  latent  heat  of  phase  transformation  in  the  solid  if  applicable 


Substituting  equations  5,  6,  7,  and  8  into  the  equation  of  energy  balance 
(2)  we  get: 


f® 


^2  FdA2-Ax 


e1(X,T1)  a5(X,T2)  wa,^)  dX  +  AH3  +  AH£  +  AH 


»  dA2  t 
+  dA2  t 


Jo  ed(X.T2)  W(X,T2)  dX 
>00 

e8(X,T2)  W(X,T2)  dX 

■  0 


(9) 
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Analysis  for  Calculation  of  Total  Emlttance  e(T)  and  the  Product  of  Total 


Emlttance  and  Absorbtance  e(T)  x  a(T)  from  Spectral  Absorbtance  ot(X,T)  Data. 

Equation  (9)  Involves  terms  containing  complicated  Integrals.  We  will 
now  present  an  analysis  to  simplify  this  equation  and  enable  us  to  calcu¬ 
late  substrate  temperature  T2> 

This  analysis  contains  two  assumptions: 

(1)  Emlttance  of  a  metal  is  independent  of  the  wave-length,  the  familiar 
grey-body  assumption  used  in  the  derivation  of  Stefan-Boltzman  Law  of 
Radiation. 

(2)  The  spectral  absorbtance  of  the  deposit  £d(X,T)  varies  very 
slowly  with  substrate  temperature,  i.e.,  for  the  purposes  and  within  the 
scope  of  this  model  it  is  essentially  independent  of  temperature. 

The  validity  of  these  assumptions  will  be  tested  by  computation 
of  the  substrate  temperature  and  comparison  with  experiment.  Making 
the  grey-body  assumption,  equation  6  which  gives  the  energy  radiated  by 
an  infinitesimal  area  of  the  deposit  dA2  becomes: 


r°° 


dA2 


eda,T2)  W(X,T2)  dX  =  ed(T2)  0  t£  dA2  t  (10) 


5  2o  A 

where  a  m  Stefan-Boltzmann  Constant  -  5.668  x  10  ergs/sec  cm  K 


Total  hemispherical  emlttance  of  the  deposit  at  substrate 
temperature  T2 


Equation  (10)  can  be  written  as 


fCO 


ed(X,T2)W(X,T2)  dX 


ed<T2> 


o  T„ 


ed(X,T2)  df  (X,T2) 


(11) 
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where 


f(A,T2> 


r°° 


c^  dX 


^  X5{exp(c2/XT2)  -  1} 


a  T. 


(12) 


The  function  f(X,T)  has  been  defined  in  reference  3. 

Since  Kirchoff's  law  states  the  spectral  emittance  is  equal  to  the 
spectral  absorptance  at  the  same  temperature,  equation  11  can  be  written 
as 


£d(T2) 


f  1 


a^X,T2)  df  (X,T2) 

0 


(13) 


Similarly, 


where , 

e.<V 


s 


(t2)  = 


aa(X,T2)  df  (X,T2) 


;0 

total  hemispherical  emittance  of  the  substrate  at  substrate 
temperature  T2« 


(14) 


as(X,T2)  ■  spectral  absorbtance  of  substrate  at  substrate  temperature. 

If  we  assume  the  variation  of  a,(X,T.)  and  a  (X,T.)  with  temperature  are 

a  l  s  4 

very  small,  i.e.  they  are  temperature  insensitive,  equations  13  and  14 
reduce  to  equations  15  and  16  respectively. 


r  1 


e„(I2)  - 

L 

ad(X) 

df  (X,T2) 

(15) 

1 

es(V  * 

0 

a8(X) 

df  (X,T2) 

(16) 

Similarly,  by  applying  the  grey-body  assumption,  QR,  the  energy 
absorbed  by  the  deposit  at  temperature  T2  from  source  radiation  at  temperature 


T1  is; 


ad(X,T2)  e1(X1,T1)  W(X,T1)  dX  dA2  t 


a(T2)  ed^  oT^  dA2  t 


(17) 


f 


1.7 


where  a(T2)  *  total  hemispherical  absorbtance  of  deposit  at  substrate 

temperature  T2 

e(T^)  ■  total  hemispherical  emittance  of  source  at  source 

temperature 

Equation  (17)  can  be^written  as  - 

|  ad(X,T2)  e1(X,T1)  W(X,T1)dX 

a(T2)  eCT^  -  -5 - * -  (lg 


Applying  Kirchoff's  law  e^(X,T)  -  a^(X,T)  and  equation  18  becomes 


a(T2)  eC^)  -  L  ad(X,T2)a1(X,T1)  Wtt.T^  dX 


Again  assuming  that  the  spectral  absorbtance  is  essentially  insensitive 
to  temperature,  the  above  equation  reduces  to 


a(T2)  eC^) 


Q  ad(X)W(X,T1)  dX 


(1 

,  0  “ 


-  j0  ad  (X)  dftt,^)  (19) 

where  ad(X)  "  spectral  absorbtance  of  deposit  of  source  material  whose 
value  is  temperature  insensitive 

and  f(X,T^)  -  the  function  defined  in  equation  12  above  at  temperature  T1 
Equations  15,  16,  and  19  can  be  used  to  simplify  the  equation  of  energy 

balance  (equation  9)  which  now  becomes 


V2  FdA2-A1t  7  T1  °d  df^.T^dX  +  Hs  +  +  Ht 

fl  0  fl  <20> 

dA2t  aT„4  ad(X)  df(X,T2)  +  dA2  t  a  T24  a8(X)  df(X,T2) 


0 


0 


Equation  20  is  much  easier  to  evaluate  than  equation  9  because  the 
integrals  in  equation  20  can  be  evaluated  using  Simpson's  Rule.  We  will 


illustrate  this  for  the  R^  term  which  is  stated  in  equation  1>  and  is 
the  first  term  on  the  right-hand  side  of  equation  20. 

The  function  f(A,T2)  in  equation  12  can  be  expanded  to 

,00 


f(A,T2) 


27thc2  dA _ 

^  A5  exp(hc/k  T2)  -  1 


(21) 


We  note  that  the  function  f(A,T2)  is  a  function  of  A  and  T2-  Its 
values  have  been  tabulated  by  Edwards  in  reference  3.  In  addition, 
the  data  for  the  spectral  absorbtance  of  Ti,  Ct ,(A)  at  room  temperature 
are  given  in  figure  2.^^  In  the  following  Table  I^we  have  tabulated 
these  values  and  applying  Simpson's  Rule  we  get 


ed(T2> 


_  .  /  A  \  1  r-  0. 8  X  2.76  A  .  A  i 

a(A)  df  =  - — - ■  0.184 

■  0 


Similarly  we  can  evaluate  the  other  integrals  in  equation  20. 


Comparison  of  Calculated  Deposit  Temperatures,  with  Measured  Values 

We  shall  calculate  the  deposit  temperature  for  the  experimental 

conditions  for  Runs  40  and  41  in  reference  1  as  shown  in  Table  II.  For 
the  experimental  geometry  used,  ,  the  configuration  factor  from 

an  infinitesimal  area  on  the  substrate  dA2  to  the  source  area  is 

0.0107.  The  thermodynamic  and  thermophysical  values  with  symbols  used 

are  listed  in  Table  III.  We  will  illustrate  the  calculation  for  Run  40. 


SPECTRAL  ABSORPTANCE 
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QR  “  dA2  FdA2-Ax  tatl 


4 


ad  (X)  dfd.T^ 


-5 


-  dA2  x  0.0107  x  3600  x  5.67  x  10  x  (2283)  x  0.095 

O 

■  53.8  x  10  x  dA2  ergs 


(23) 


Qc  Is  evaluated  by  cooling  the  metal  vapors  from  source  temperature  2283°K 
to  1735°K,  condensing  them  to  solid  at  1735°K  and  cooling  the  solid  to  the 
substrate  temperature  T2  .  Now, 

Qc  "  Hs  +  H*  +  Ht 


H  - 

s 


,  2283 


1735 


,1735 


dT  + 


C  dT  /  p  6  dA2 


(24) 


-  [548  x  5.1  +  (1735  -  T2)]  x  106  x  4.6  x  0.094  x  dA£  ergs 
where  p  -  density  -4.6  g/cm-*  for  Ti  and  6  -  thickness  of  Ti  deposit 
near  center  ■  0.094  cm. 


7735°K)  *  L  p  6dA2  “  9.72  x  10^  x  4.6  x  0.094  x  dA2  ergs  (25) 

■  energy  associated  with  the  latent  heat  of  transformation  of  a  to  8  at 
1155°K 

*  8.86  x  108  x  4.6  x  0.094  x  dA2  -  3.83  x  108  ergs  (26) 

Rd  "  dA2  °  T2^Ed^T2^ 

■  dA2  x  3600  x  5.67  x  10  5  x  x  0.184 

-2  4 

*  3.76  x  10  x  T2  x  dA2  ergs  (27) 

*s  ■  “2  '  0  T24  es(T2) 


Substituting  the  value  of  eg(T2)  for  stainless  steel  as  0.22 
R  -  dA„  x  3600  x  5.67  x  10"5  x  T*  x  0.22 

8  Z  ^ 

-  4.49  x  10”2  x  x  dA2  ergs  (28) 

Substituting  equations  23,  24,  25,  26,  27,  28  into  equation  20  we 

4 

get  a  biquadratic  equation  containing  a  T2  and  T2  term. 
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Solving  for  T£  we  get  T£  ■  895  °K  for  Run  40 
Similarly  we  calculate  T2  ■  913  °K  for  Run  41. 

Table  IV  compares  the  calculated  the  experimental  values  of  the  deposit 
temperatures  and  the  percent  difference  in  values.  The  agreement  is 
excellent,  the  difference  being  only  2°K  or  .02%  for  Run  40  and  26°K 
or  less  than  3%  for  Run  41.  Thus  we  may  conclude  that  the  model 
postulated  and  the  simplifying  assumptions  made  to  calculate  the  substrate 
temperature  are  reasonable  and  well  within  the  uncertainty  bounds  of  the 
input  data  (precision  of  measured  source  and  substrate  temperatures, 
precision  of  thermodynamic  data  used  etc.) 

It  would  be  instructive  to  compute  the  error  in  the  deposit  tempera¬ 
ture  for  deviations  of  ed(T2) .  eg(T2)  and  ad^T2^  x  E1^T1^  from  the  values 
used  in  the  original  computation  of  Tj.  This  is  done  in  two  ways.  First, 
it  is  assumed  that  all  three  parameters  ed(T2),  eg(T2)  and  ad(T2)  x  e^Tj) 
are  off  by  +  10Z,  +  20Z,  -  10Z  and  -  20Z.  The  results  are  shown  in  Table  V. 
We  note  that  the  largest  deviation  of  the  computed  T2  value  Is  only  6Z. 
Secondly,  it  was  assumed  that  only  e^Tj),  the  total  hemispherical  emittance 
of  liquid  titanium  is  off  by  +  10Z,  +  20Z,  -  10Z,  -  20Z.  The  computed 
substrate  temperature  values  are  only  slightly  effected. 

Calculated  Temperature  Distribution  of  Deposit 

To  compute  the  temperature  distribution  from  center  to  edge  of  the 
deposit,  we  need  to  know  the  appropriate  configuration  factors  and 
thickness  of  the  deposit  at  various  points  and  substitute  them  into  the 
equation  of  energy  balance.  Using  the  measured  thickness  of  the  deposit 
in  Run  40  at  various  points,  the  deposit  temperature  distribution  is 
computed  and  tabulated  in  Table  VII. 
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Summary  and  Conclusions 

A  model  has  been  proposed  for  calculating  the  temperature  of  the 
substrate  during  vacuum  deposition.  It  Is  based  on  an  energy  balance  In 
the  substrate-deposit  combination.  The  energy  Input  Is  from  thermal 
radiation  from  the  evaporation  source  and  from  the  condensation  of  metal 
vapors.  The  energy  loss  is  by  thermal  radiation  from  the  substrate 
and  the  deposit.  Two  assumptions  were  made  to  simplify  the  equation 
of  energy  balance  to  permit  computation  of  substrate  temperature  and 
comparlsion  with  experiments.  They  were  -(1)  The  "Grey  Body"  assumption, 
l.e.  spectral  emlttance  Is  Independent  of  wave-length  and  -(2)  spectral 
emlttance  was  relatively  Insensitive  to  temperature.  Computation  of 
substrate  temperature  agreed  very  well  with  the  experimentally  measured 
temperature  for  titanium  evaporation.  Thus  the  model  and  simplifying 
assumptions  appear  to  be  valid. 

It  Is  proposed  that  with  this  model  one  can  set  up  evaporation 
experiments  whereby  measuring  the  substrate  temperature  and  other 
experimental  parameters,  the  total  hemispherical  emlttance  of  liquid 
metals  can  be  obtained.  Such  data  Is  very  sparse  for  liquid  metals. 
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Appendix  1;  Heat  transfer  due  to  conduction 

We  have  neglected  the  mode  of  energy  transfer  due  to  conduction  in 
the  equation  of  energy  balance.  The  following  calculation  will  show 
that  this  term  is  small  compared  to  that  due  to  radiation. 


For  experiment  Run 

No. 

40,  we  have  previously  calculated  - 

-2  4 

3.76  x  10  x  T2 

dA2 

ergs 

(27) 

4.49  x  10  2  x  T^ 

dA2 

ergs 

(28) 

We  consider  an  area  of  5  cm.  radius  at  the  center  of  the  deposit  and 
assume  that  the  substrare  and  deposit  radiant  at  an  average  temperature 
of  1/2(895  +  877)  ■  886°K  where  895°K  is  the  center  temperature  and  8776K 
is  the  edge  temperature  of  this  5  cm  circular  area.  Thus,  the  total 
energy  radiated  is 

Rd  +  Rs  -  (8.25  x  lO-2)  x  (886)4  x  n  x  25  -  4  x  1012  ergs  (1') 

We  define  the  following  terms 

■  energy  transferred  by  conduction  from  the  center  to  points  5  cm 
away  in  the  deposit 

Eg  «  energy  transferred  by  conduction  from  the  center  to  points  5  cm 
away  In  the  substrate. 


Now, 

E  .  -  -  K .  A .  t 
d  d  dr  d 


where 


(2’) 


K.  *  the  thermal  conductivity  of  Ti 


dT 

dr 


■  1.72  x  10  ergs/sec  cm  °K 

■  the  temperature  gradient  from  center  to  points  5  cm  away 

-  - 5 -  ”  3.6  k"/cm 

■  the  circumferential  area  of  the  deposit 

0.094  +  0.086 

■  271  x  5  x  - 2 - 

-  2.82  cm2 
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t  *  duration  of  the  experiment 
■  3600  sec. 

Substituting  all  the  values  into  equation  (2')  we  get 

E.  *  -  1.72  x  106  x  3.6  x  2.82  x  3600 
a 

9 

=  -  62.9  x  10  ergs 
Similarly! 


E  *  - 


s 

where 


K 


8 


t 


Ka  ■  the  thermal  conductivity 
*  2.84  x  10^  ergs/sec  cm  "K 


O') 

(4’) 


A  =  the  circumferential  area  of  stainless  steel 
8 

-  2ir  x  5  x  0.0127 

n  /  2 

-  0.4  cm 

Substituting  all  values  into  equation  (4 * > ,  we  get 

E  -  -  2.44  x  10^  x  3.6  x  0.4  x  3600 
s 

-  -  11.7  x  10^  ergs  (5') 

Thus,  the  total  energy  transfer  from  the  center  to  points  5  cm  away  due 
to  conduction  is  E^  +  Eg  *-7.46  x  10^  ergs. 

Since  it  is  less  than  2X  of  the  energy  due  to  radiation  of  the 
substrate  and  deposit,  we  are  justified  in  neglecting  it. 


TABLE  I 


Application  of  Simpson's  Rule  to  Evaluation  of  £^2) 


f 

he 

kXT 

XT 

X  at  893°K 

ad(X) 

Weighting 

Factor 

Product 

.1 

6.55 

2200 

2.46 

0.31 

1 

0.31 

.3 

4.61 

3120 

3.50 

0.26 

4 

1.04 

.5 

3.50 

4110 

4.61 

0.23 

2 

0.46 

.7 

2.58 

5580 

6.25 

0.20 

4 

0.80 

.9 

1.53 

9410 

10.54 

0.18 

1 

0.15 

+2.76 


Hence 

f1 

€(T2)  - 

o(X)  df 

0 

0.8  x  2.76 

12 

0.184 


Experimental  Condition  for  Runs  40  and  41 


Run  No. 

40 

41 

Melt  temperature 

2010°C 

2010°C 

Duration  of  experiment 

60  min. 

60  min. 

Size  of  evaporant  stock 
(diameter  and  height) 

3x4  in2 

3x5  in2 

Collector  size 

12x12  in2 

12x12  in2 

Collector  temperature 

620°C 

640°C 

Distance  of  collector 
from  source 

14.5  in 

14.5  in 
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TABLE  111 

Thermodynamic  and  Thermophysical  data  used  In  Calculations 


Symbol 

Measuring 

Value 

C 

Average  heat  capacity  at  constant 

5.1  x 

106 

ergs/gm°K 

P1 

pressure  for  Tl  vapor  between 

2283  and  1735°K 

C 

Average  heat  capacity  at  constant 

6.6  x 

106 

ergs/gm°K 

p2 

pressure  for  solid  Ti  between 

1735  and  733°K 

L 

Latent  heat  of  sublimation  of 

9.72  x 

10^  ergs/gm°K 

s 

Ti  at  1735°K 

Lt 

Latent  heat  of  a-8  phase  trans¬ 

8.86  x 

O 

00 

ergs/gm°K 

L 

formation  of  Ti  at  1153°K 

Reference 

6 

6 

6 

6 


es(T2)  Total  hemispherical  emlttance  0.22 

of  stainless  steel  at  tempera¬ 
ture  800-1000°K 


Calculated  In 
this  paper 


ed(T2)  Total  hemispherical  emlttance 
of  solid  Ti  at  temperature  T2 


0.184 


Calculated  In 
this  paper 


ad(T2)  e^(T^)  Product  of  total  hemispherical  0.095 

absorptance  at  substrate  tempera¬ 
ture  T2  and  the  total  hemispherical 
emlttance  at  source  temperature  T^ 


Calculated  in 
this  paper 
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TABLE  IV 

Comparison  of  Calculated  Experimental  Values  of  Substrate 
Temperature  for  Titanium  Deposition 


Run  40 


Calculated 

895 


Substrate  Temperature  °K 


Measured 

893 


Difference 

2 


913 


26 


%  Difference 

0.02 


Run  41 


887 


2.6 


1.19 


TABLE  V 

Calculated  Substrate  Temperature  for  Various  %  Deviations  In  Values  of 

EjttO.  W.  0tj(T2)  X  E1(T1)  for  Runs  40  and  41  -  °K» 


Experimentally  Calculated  Values  for  %  Deviations  of 


Measured 

0% 

-10% 

-20% 

+10% 

+20% 

Run  40 

893 

895 

915 

938 

878 

865 

Run  41 

913 

887 

906 

928 

871 

856 

TABLE  VI 

Calculated  SubBtrate  Temperatures  T»  for  Various  %  Deviations  In  Value  of 
e^(T^)  -  the  Total  Hemispherical  Emlttance  of  Liquid  T1  source  -  °K 


Experimentally 


Calculated  Values  for  %  Deviations  of 


Measured 

0% 

-10% 

-20% 

+10% 

+20% 

Run  40 

893 

895 

893 

890 

898 

900 

Run  41 

913 

887 

884 

882 

888 

891 
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TABLE  VII 

Temperature  Distribution  of  Deposit 


Distance  away  from  center 
0.001  cm 
5  cm 
10  cm 
15  cm 


Experimental  Value 
of  the  thickness 

0.094  cm 

0.086  cm 

0.066  cm 

0.043  cm 


Temperature 

895°K 

377°K 

825°K 

748°K 
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SUPPLEMENT-2 

Temperature  and  Thickness  Distribution  on  the  Substrate  During 
High  Rate  Physical  Vapor  Deposition  of  Materials 

I.  INTRODUCTION 

High  rate  physical  vapor  deposition  (HRPVD)  is  a  relatively  new 
technique  for  synthesis  of  metals,  alloys,  ceramics,  dispersion  strengthened 
alloys,  etc/1’2,3*  Deposition  rates  of  .001  inch  (250,000  A)  per  minute 
or  greater  are  obtained.  The  products  which  are  self-supported  structures 
or  thick  coatings  have  practical  applications  in  advanced  engineering 
structures  such  as  honeycomb  core  material,  coatings  for  Jet  engine  turbine 
blades,  etc.  Essentially,  the  technique  consists  of  two  processes; 
evaporation  and  condensation  .  Metals,  alloys  and  ceramics  contained  in 
water  cooled  copper  crucibles  are  evaporated  in  high  vacuum  by  using  high 
power  electron  beams.  The  vapors  are  collected  on  heated  metallic  sub¬ 
strates.  The  deposition  temperature  is  a  very  important  process  parameter. 

It  is  desirable  to  obtain  a  material  of  full  density  and  small  grain  size 

(4) 

for  high  strength  and  toughness.  Bunshah  and  Juntz  have  made  the  following 
observations.  At  high  deposition  temperatures,  the  material  is  fully  dense  but 
has  a  large  grain  size.  As  deposition  temperature  decreases  the  grain  size 

becomes  smaller.  At  very  low  deposition  temperatures,  the  material  is 
less  than  full  density.  Hence  it  becomes  important  to  deposit  at  an 
optimum  temperature  to  obtain  fine  grain-size  fully  dense  deposits. 

In  the  simple  deposition  geometry  with  a  flat  substrate  placed  per¬ 
pendicular  to  the  axis  of  the  vapor  source,  see  Fig.  1,  the  thickness  of 
the  deposit  varies  from  the  center  to  the  edges.  The  thickness  gradient 
and  the  rate  of  deposition  decrease  as  the  source-substrate  distance  increases. 
Therefore,  there  is  an  optimum  source-substrate  distance  to  obtain  the 
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Figure  1 


Schematic  of  High  Rate  Physical  Vapor  Process 
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maximum  yield  of  uniform  thickness  (/V  5%)  deposit  consistent  with  a  high 
deposition  rate. 

The  experimental  variables  in  the  high  rate  physical  vapor  deposition 
process  are  as  follows: 

1.  Nature  of  the  evaporating  species  (metal,  alloy,  compound)  and  the 
vapor  pressure  of  the  component  species. 

2.  Temperature  of  the  source  T^  and  the  evaporation  rate. 

3.  Thermodynamic  properties  of  condensate  (latent  heat  of  condensation, 
spaciflc  heat,  vapor  pressure). 

4.  Heat  of  reaction  In  the  case  of  reactive  evaporation,  e.g.  metal  vapor 
plus  oxygen  readtlng  to  give  metal  oxide.  The  reaction  may  take  place 

on  the  substrate  or  in  the  gas  phase. 

5.  Source  diameter. 

6.  Source  to  substrate  distance,  H. 

7.  Substrate  material. 

8.  Substrate  thickness. 

9.  Size  of  the  substrate* 

10.  The  total  hemispherical  emittance  of  the  evaporant  as  a  function  of 
evaporation  temperature,  E^(T). 

11.  The  total  hemispherical  emittance  of  the  deposit  as  a  function  of 
condensation  temperature,  e^CT) . 

12.  The  total  hemispherical  emittance  of  the  substrate  as  a  function  of 

condensation  temperature,  E  (T) . 

s 

13.  The  absorptance  of  the  deposit  as  a  function  of  condensation  tempera¬ 
ture,  Oj(T) . 

14.  The  specific  heat  of  the  substrate  material  as  a  function  of 
temperature . 

15.  The  specific  heat  of  the  condensate  material  as  a  function  of 
temperature. 


II.  OBJECTIVE  OF  THIS  INVESTIGATION 
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With  the  large  number  of  process  variables  it  is  important  to  be 
able  to  focus  attention  on  the  range  of  these  variables  to  be  explored 
experimentally  with  the  objective  of  obtaining  the  desired  structure  in 
the  deposit.  To  this  end  it  is  the  objective  of  this  investigation  to 
illustrate  a  method  for  calculating  the  following: 

1.  Deposition  temperature. 

2.  Deposition  temperature  distribution  along  the  substrate. 

3.  Deposit  thickness  and  its  variation  along  the  substrate. 

The  deposition  of  yttrium  and  yttria  (YjO^)  hy  evaporation  from  a 
liquid  yttrium  source  will  be  used  to  illustrate  the  method.  ^2^3  *8 
formed  by  reactive  evaporation  of  Y  vapors  in  the  gas  phase  or  on  the 
substrate  with  introduced  into  the  gas  phase.  The  results  of  the 
calculations  will  yield  the  specific  range  of  variables  to  be  explored 
experimentally. 

Radiation  from  the  source,  energy  liberated  by  the  condensation  of 
metal  vapors  and  the  heat  of  reaction  for  reactive  evaporation  are  the 
sources  heating  the  substrate.  Their  relative  importance  for  various 
experimental  conditions  will  be  evaluated. 

III.  CALCULATIONS 

Calculations  will  be  made  for  the  following  range  of  experimental 
variables. 

1.  Temperature  of  the  liquid  Y  source  2,200,  2,300,  2,400°K  (M.P.  of 
Yttrium  is  1768°K) • 

2.  Source  diameter  -  1.0  and  1.5  inches. 

3.  Source  to  substrate  distance  -  10,  15,  25,  50  cm. 


4. 


Substrate  materials  -  copper,  stainless  steel. 
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5.  Substrate  thickness  -  0.001,  0.005,  0.01,  0.05  inch, 

6.  Substrate  size  (diameter)  -  8,  16,  24  cm. 

A.  Vapor  Pressure  and  Evaporation  Rate 

Vapor  pressure  data  are  given  in  the  following  form 

Log1()  P  -  A  T_1  +  B  lo8lOT  +  CT  +  DT2  +  E 

where 

P  ■  vapor  pressure  of  the  metal,  torr 
T  ■  temperature  of  the  metal,  °K 

A,B,C,D,  and  E  -  constants,  characteristic  of  the  metal 
Evaporation  Rate 

The  evaporation  rate  of  a  metal  is  governed  by  Langmuir  equation 

G  -  K  P  (M/T)1/2 

where 

2 

G  ■  weight  evaporated,  g/cm  -sec. 

P  ■  vapor  pressure,  torr. 

T  ■  temperature  of  the  metal,  eK. 

M  -  molecular  weight  of  the  gas  phase 

K  ■  constant  *  0.05833 

The  vapor  pressure  and  evaporation  rate  of  yttrium  at  various 
temperatures  are  calculated  and  tabulated  in  Appendix-1.  Evaporation  rate 
of  yttrium  with  temperature  is  plotted  in  Fig.  2. 

B.  Temperature  and  Temperature  Distribution  of  the  Substrate 

The  experimental  set-up  for  high  rate  plysical  vapor  deposition 
processes  (HRPVD)  is  shown  in  Figure  1,  The  metal  is  contained  in  a 
water-cooled  copper  crucible  heated  by  an  electron  beam  source.  The 
metal  vapors  are  condensed  on  a  metallic  substrate  located  above  the 
source . 


EVAPORATION  RATE  (G/CM2-  SEC) 


Figure  2a.  Evaporation  Rate  of  Yttrium  vs  Temperature 
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The  model  for  calculation  of  the  substrate  temperature  Is  based  on  the 

energy  balance  between  the  energy  input  and  output  of  the  substrate. 

The  energy  Input  to  the  substrate  consists  of  three  terms: 

Q_  -  the  radiant  energy  from  the  source  at  temperature  T 
K  '  1 

Qc  -  energy  supplied  by  the  condensation  of  metal  vapors 

-  the  heat  of  reaction  for  reactive  evaporation  i.e.,  if  a  reaction 
occurs  between  the  condensing  vapors  and  a  gas  phase  on  the  substrate 
or  between  two  metal  vapors  condensing  on  the  substrate.  For  simple 
condensation  of  metal  or  alloy  vapors  this  term  is  zero. 

The  energy  output  from  the  substrate  consists  of  the  following  terms*. 

Rj  -  energy  radiated  from  the  deposit  at  temperature  of  the 

substrate 

R  -  energy  radiated  from  the  back  side  of  the  substrate  at  temperature 
s 

Tj  i.e.,  the  side  away  from  the  source 

H  -  energy  necessary  to  heat  the  substrate  from  the  initial  temperature 
s 

T^  to  the  final  temperature  T^  in  any  interval  of  time 

Hd  -  energy  necessary  in  the  nth  interval  to  heat  material  deposited  in 
the  (n-1)  intervals. 

We  can  now  write  the  equation  of  energy  balance  as: 

\  +  «C  +  V  ■  Rd  +  R8  +  Hs  +  Hd  (1) 

Initially  during  the  process,  the  energy  input  will  exceed  the  heat 

output  and  the  temperature  of  the  substrate-deposit  combination  will  rise; 

thus  causing  an  increased  energy  loss  from  the  substrate.  The  starting 

temperature  of  the  substrate  is  assumed  to  be  300°K.  Next,  time  is  divided 

into  intervals  of  10  seconds  and  using  equation  1  above,  the  temperature 

rise  in  the  substrate-deposit  combination  is  computed  for  each  interval. 

The  final  temperature  of  each  interval  is  taken  as  the  starting  temperature 

of  the  following  one.  Since  the  energy  input  to  the  substrate  deposit 

4 

combination  is  constant  and  the  output  is  proportional  to  T^  (T£  being 
the  substrate  temperature) ,  a  steady  state  will  be  reached.  Under  steady 


2.6 


state  conditions,  H  and  in  equation  1  are  zero  and  the  equation  of 
energy  balance  becomes 

«R+«c  +  V  Rd  +  Rs  <2) 

The  validity  of  this  model  has  been  confirmed  by  Chow  and  Bunshah^  by 
comparing  the  computed  and  measured  substrate  temperatures  for  titanium 
evaporation.  In  further  calculations,  the  following  assumptions  were  made: 

1.  The  surface  of  the  evaporating  source  is  flat  and  at  a  uniform 
temperature. 

2.  Heat  loss  by  conduction  on  the  substrate  is  negligible,  as  shown  by 
Chow  and  Bun  shah. 

3.  The  surface  of  the  deposit  behaves  like  a  gray  body  and  hence 
spectral  absorptance  is  equal  to  the  spectral  emittance.  This  assumption 
is  often  used  in  radiant  heat  transfer  calculations^.  Its  validity  has 
been  tested  by  Chow  and  Bunshah^  for  the  deposition  of  titanium. 

4.  The  substrate  has  black  body  surroundings  at  a  constant  temperature 
of  300°K. 


5.  The  deposit  has  theoretical  density. 

Let  us  now  examine  in  detail  each  of  the  terms  in  Equation  (1)'. 

Qr  -  Energy  gain  due  to  radiation  from  the  source. 

Heat  absorbed  by  an  infinitesimal  area  dA^  on  the  substrate  due  to 
radiation  from  the  source  of  area  and  using  the  reciprocity  relationship 
shown  later  in  eq.  9,  we  get: 

4  <3> 

qr  -  a  dA2  El(T>  ad(T)  I2 

where 


VT> 


4  2 

Stef an -Boltzmann  Constant,  cal/°K  -sec. -cm. 

Total  hemispherical  emittance  of  the  source 
Total  hemispherical  emittance  of  the  deposit 
Total  hemispherical  absorptance  of  the  condensate 
Temperature  of  the  source,  °K. 
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FdA  -A  =  Conf iguration  factor  for  radiation  going  from  area  dA2  to  area  A^. 

It  is  assumed  that  the  absorbtance  of  the  deposit  has  to  be  taken 
as  that  of  the  substrate  until  50  atomic  or  molecular  layers  are  deposited. 

The  times  to  achieve  this  condition  for  the  experimental  conditions  under 
investigation  have  been  calculated  and  given  in  Appendix-2.  The  times  are 
too  small  compared  to  the  duration  of  the  experiments  and  hence  are 
neglected. 

Qc~EnerRy  gain  due  to  condensation  of  vapors 

Energy  given  to  an  infinitesimal  area  dA2  on  the  substrate  due  to 
condensation  of  vapors  is: 

Q  -  number  of  moles  condensing  x  heat  of  condensation  per  mole. 

c 

The  number  of  moles  condensing  per  second 


=  M  dA2  FdA2-A1 
where  M  ■  molecular  weight 

2 

G  -  evaporation  rate,  g/cm  sec 
Heat  of  condensation  per  mole 


■  (H  -  K_no)  +  heat  of  sublimation  at  298°K 

298  vapor 

-  heat  required  to  heat  Y  from  298°K  to  the  substrate  temperature  T2 

where 


(H„  -  H_no)  ■  enthalphy  change  in  vapor  from  T,  to  298CK. 

Tj  298  vapor  1  °  1 

Let 


F  ■  (H_  -  Hon0)  +  heat  of  sublimation  at  298°K. 

T^  298  vapor 

The  values  of  F  are  computed  as  a  function  of  source  temperature  and  given 
in  Appendix  3. 

Heat  required  to  heat  one  mole  of  yttrium  from  298°K  to  the  substrate 
temperature  ?2 

=  molar  specific  heat  x  (T2  -  298) 
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_<*  0.29  x  10  .  _ 

=  (5.59  +  1.9  x  10  x  T2  +  2  >  (T2  "  298) 

«=  5.0238  x  T2  +  1.9  x  10"3  x  T22  +  -,29  x  10 


86.42  x  10' 


-  1665.82 


Q. 


G  dA2FdA2-A1 


M 


{  F  -  5.0238  x  T„ 


.  ^  . „-3  „  2  0.29  x  105  .  86.42  x  105  .  nn1  ,  , 

-1.9x10  x  T2  -  - - -  + - 2 - +  1665 . 82}  cal/ sec . 


-  Heat  gain  due  to  Reaction 

The  model  assumed  for  the  deposition  of  yttria  is  that  yttrium  vapors 
condense  on  the  substrate  and  then  react  with  gaseous  oxygen  to  form  solid 
yttria: 


vapcr 


solid 


3 

2Y  +  ~  0 

(solid)  2  U2 


- ►  l2°3 

(gas)  J (solid) 


Since,  heat  of  reaction  is  a  thermodynamic  state  property,  even  if  metal 

vapors  react  with  oxygen  gas  in  the  vapor  phase  to  form  oxide  and  then 

condense  to  solid  yttria,  the  heat  of  reaction  will  be  the  same. 

Heat  given  to  an  infinitesimal  area  dA  on  the  substrate  due  to 

2 

the  formation  of  yttria  : 


<Rx 


number  of  moles  of  yttria  formed  x  heat  generated  per  mole. 


The  number  of  moles  of  yttria  being  formed  per  second 

°  r«^i 


2M 


Heat  generated  per  mole  is  equal  to 


ah298ok  + 


r  t. 


ACp  dT 


(4) 


298 
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where 


AH2gg  =  Enthalpy  of  reaction  at  298°K,  cal/mole. 
AC  =  C 


-  2  c  -  r 

P  PV  0  Py(s)  2  P 


Y2°3 


2  (g) 


=  7.68  -  4.1  x  10'3  x  T.  -  10‘ 

1  „  2 


AC  =  7.68  xT-  2.05  x  lo"3  x  T  2  +  ±L7?-£..3t)5  -  3706.60 
298  P  1  T2 


G  dA2  FdA  -A 

QRx  =  - 2M~  ~  |^459‘  106-60  _  7<68  x  T2  +  2*05  x  10"3  x  T22  "  10^ 

cal/sec.  (5) 


R  -  Heat  Loss  Due  to  Radiation  from  Substrate 
-s  — - — 


The  heat  loss  due  to  radiation  from  the  substrate  to  the  surroundings, 


from  an  area  dA2> 


where 


R  =  a  e  (T)  (T  4  -  3004)  dA„  cal/sec. 

S  8  L  L 


(6a) 


e  (T)  -  total  hemispherical  emittance  of  the  substrate 

3 

Rj  -  Heat  Loss  due  to  Radiation  from  the  Deposit 

The  heat  loss  due  to  radiation  from  the  deposit  to  the  surroundings, 
from  an  infinitesimal  area  dA2 , 


Rd  "  aed^T^  ^2  ~  300^  flA2  cal/sec* 


(Cb) 


Hfl  -  Heat  required  to  heat  the  substrate  from  the  starting  temperature 
(Tj)  of  an  interval  to  the  final  substrate  temperature  (T7|)  of  the  interval. 

The  heat  required  to  heat  an  element  of  substrate  of  area  dA2  from 
the  starting  temperature  of  the  interval  (T^)  to  the  final  temperature 
of  the  interval  (T^) , 

Hs  ■  mass  x  specific  heat  x  (T^-T^) 

*  thickness  x  density  x  dA2  x  specific  heat  x  (T^-T^)  (7) 
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th 


Hj  -  Heat  required  In  the  n  Interval  to  heat  the  condensate  deposited  in 
the  prior  (n-1)  Intervals  from  the  starting  (T^)  to  the  final  (T^)  tempera- 
ture  of  the  interval 


For  the  deposition  of  yttrium, 

Hj  =  number  of  moles  of  Y  deposited  in  (n-1)  intervals  in  an  area  dA2 
x  molar  specific  heat  of  Y  x  (T^-T^) 

=  £  dA_  F..  .  (n-1)  [5.59  x  T.  +  1.9  x  10_3  x  T  2 

M  2  “  4  (8a) 

.  0.29  x  10*  x  T_ 

+  -°-‘-2|— -  5.59  x  I-  -  1.9  x  10~3  x  T  x  T, - x - -j 

i.  a  ^  3  ^  ^ 


For  the  deposition  of  yttria, 

=  number  of  moles  of  yttria  deposited  in  (n-1)  intervals  in  an  area 


dA. 


x  molar  specific  heat  of  yttria  x  (T^  -  T^) 


G  d/2  FdA  -A  (n_1)  I29-60  x  T4  +  1*20  x  10  3  x  T42 


2M 

4.78  x  10 


-  29.60  x  T3  -  1.20  x  10“3  x  T3  x  T4 


4.78  x  10J  x  T_ 

+  - o - h 


(8b) 


Computation  of  the  Substrate  Temperature  Tp 

Substituting,  the  above  equations  3, 4, 5, 6, 7  ,8  representing  the  various 
terms  into  the  equation  of  energy  balance  (eq.  1),  it  assumes  the  form: 

A  A 

T2  +  A5  X  T2  "  IT  +  A7  X  T2  +  tT  “  A10  ° 

Art  * 


where,  A^  A^,  A7>  Ag  and  A1Q  are  material  constants. 


,10 


An  equation  of  this  form  is  solved  for  T2  by  Newton-Ralphson  method  which 
is  as  follows! 


f  “  T2  +  A5  X  T2  "  T  2  +  a7  x  t2  +  t2  _  A10 

2  2  x  A,  A 

f'  *  4xT2  +  2  x  Aj  x  T2  +  - j-  +  A^  -  — 2 

T2  T2 
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Initially,  T?  is  given  an  arbitrary  value,  say  T  '  .  Then  a  new  value  of  T,  is 

found,  e.g.,  T2  new  =  — (f / f * ) .  This  new  value  is  substituted  in 

place  of  T^  and  another  new  value  is  /ound.  This  process  is  repeated 

until  T  new  -  T'  is  within  the  limits  of  accuracy  desired  (1.0%  is 
2  2 

sufficiently  accurate  for  cur  calculations). 

A  computer  program  has  been  developed  to  compute  the  temperature 
on  the  substrate  at  various  points,  as  a  function  of  experimental  variables. 
The  temperature  variation  is  plotted  in  figures  3  &  4  and  tabulated  in  tables 
I  and  II. 

The  data  used  in  these  calculations  is  given  in  Table  XXII  along  with 
the  symbols  and  sources  of  data.  Since  no  values  were  available  for  the 
total  hemispherical  emittance  of  Y,  the  spectral  value  at  \  =  0.65M  is  used 
instead.  Similarly  for  the  total  hemispherical  absorbtance  of  Y^,  the 
spectral  value  at  A  =  0.65y  for  oxidised  yttrium  is  used.  Changes  in  the 
computed  substrate  temperature  T2  for  deviations  from  the  assumed  values 
of  emittance  and  absorbtance  are  calculated  and  discussed  below. 


C.  Thickness  Distribution  on  the  Substrate 

The  mass  of  yttrium  reaching  an  infinitesimal  area  dA2  on  the 
substrate,  from  the  source  is  equal  to; 


G  Ai  dFA 

1  A^-dA2 

where 

2 

A  •,  =  surface  area  of  the  source,  cm 
Using  the  reciprocity  relation  of  configuration  factors, 


G  A1  dFArdA2  “  G  dA2  FdA2-Al 


where 


(9) 


dF 


Ax-dA2 


dA2"Al 


configuration  factor  to  dA2  from  the  source  of  surface 
area  A. 


SUBSTRATE  MATERIAL-COPPER 
TEMPERATURE  CF  THE  SOURCE  (°K) 


Figure  3.  Temperature  Variation  for  the  Deposition  of  Yttrium 
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where  Z  =  1  +  (1+Y2)  xl 
Y  =  b/H 
x  =  H/6 

and  6,  b,  and  H  are  shown  In  figure  1. 

The  thickness  at  any  point  on  the  substrate  Is  calculated  by  using  the  mass 
balance  equation 


where  T 
P 

T 


G  dA2  FdA2-A1  =  T 


dA2  p 


thickness  of  the  condensate  at  any  point,  cm 
density  of  the  deposit,  g/cnf* 


dA2-Ax 


cm. 


The  calculated  thickness  values  are  shown  in  figure  5  and  Table  III. 

IV.  DISCUSSION  OF  RESULTS 


(10) 


The  results  of  the  computations  will  be  discussed  under  four  principal 
headings: 

A.  Effect  of  Experimental  Variables  on  Temperature  and  Thickness  Distribution 
on  the  Substrate. 

B.  Time  to  Reach  Steady  State  Temperature  as  a  Function  of  Substrate 
Thickness. 

C.  Errors  in  Substrate  Temperature  Calculations  due  to  Deviations  in  the 
Assumed  Emittance  Values  for  Source  and/or  Substrate  Materials  and  of  the 
Surroundings. 

D.  Relative  Contribution  of  Various  Energy  Scurces  in  Heating  the  Substrate. 
The  discussion  is  as  follows: 

A.  Effect  of  Experimental  Variables  on  the  Temperature  and  Thickness 
Distribution  on  the  Substrate. 

1.  Source  Temperature  -  Referring  to  tables  I,  II,  &  III  as  the 
temperature  of  the  source  increases,  the  evaporation  and  the  deposition  rates 
of  the  metal  increase.  Thus,  the  deposition  rate  at  2300°K  is  2.5  times 
and  that  at  2400°K  is  5.5  times  the  deposition  rate  at  2200°K.  The  tempera¬ 
ture  on  the  substrate  increases  with  increasing  source  temperature.  The 
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corresponding  temperatures  for  the  deposition  of  yttria  are  higher  than 
those  for  the  deposition  of  yttrium.  The  percent  thickness  and  temperature 
variation  from  the  center  to  the  edges  of  the  deposit  is  independent  of 
the  source  temperature.  The  percent  temperature  variation  is  only  about 
1%  higher  for  the  deposition  of  yttria,  than  for  yttrium. 

2.  Source  to  substrate  distance  -  As  the  distance  between  the 
source  and  the  substrate  is  increased,  the  deposition  rate  and  the  temperature 
at  each  point  of  the  substrate  decrease  considerably.  But,  the  percent 
thickness  and  temperature  variation  also  decreases.  In  other  words,  as 

the  source-substrate  distance  is  increased,  the  deposit  becomes  more  and 
more  uniform  in  temperature  and  thickness  and  hence  in  structure  and  mechanical 
properties.  The  temperature  variation  along  the  deposit  for  both  yttria 
and  yttrium  is  about  same,  as  seen  in  tables  1,  II,  &  III.  The  corresponding 
substrate  temperatures  for  yttria  are  higher  than  those  for  yttrium. 

3.  Diameter  of  source  -  The  thickness  and  temperature  at  each  point 
on  the  deposit  increase  with  increasing  diameter  of  the  source.  Thus,  from 
table  III,  the  deposition  rate  with  1.5"  source  is  twice  that  with  a 

1.0"  diameter  source.  But  the  percent  thickness  and  temperature  variation 
on  the  substrati*  is  independent  of  the  source  diameter. 

4.  Substrate  size  -  As  the  diameter  of  the  substrate  increases, 
the  thickness  and  temperature  difference  between  the  center  and  the  edges 
increase.  Thus,  the  thickness  variations  from  the  center  to  the  edge  are 
25,  60,  and  80%  for  8,  16,  and  24  cm.  diameter  substrates  respectively,  for 

a  source  to  substrate  distanceyH,of  10  cm.  Of  course,  this  variation  becomes 
less  with  increasing  H.  Similarly,  as  can  be  seen  from  tables  I  &  II  the 
percent  temperature  difference  between  the  center  and  the  edges  increases  with 
the  substrate  size  for  the  deposition  of  both  yttrium  and  yttria. 
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5.  Substrate  material  -  The  temperature  at  each  point  with  a 
copper  substrate  is  higher  than  that  with  stainless  steel  substrate.  This 
is  due  to  the  lower  emittance  of  copper  as  compared  to  stainless  steel. 

The  percent  temperature  difference  for  the  two  substrate  materials  increases 
with  increasing  substrate  temperature  and  varies  from  4  to  10%,  for  the 
deposition  of  yttrium.  It  would  be  higher  for  yttria  deposition  because 
of  higher  substrate  temperatures.  The  results  are  given  in  Tbble  IV. 

B.  Time  to  Reach  Steady  State  Temperature  as  a  Function  of  Substrate  Thickness 

The  time  to  reach  steady  state  temperature  increases  markedly  with 
increasing  substrate  thickness,  as  can  be  seen  in  tables  V  &  VI.  The  time 
decreases  considerably  with  increasing  source  temperature,  but  increases 
with  increases  H.  The  corresponding  times  for  yttria  are  much  smaller 
than  for  yttrium,  especially  at  greater  su'  _ trate  thicknesses. 

C.  Errors  in  Substrate  Temperature  Calculations  due  to  Deviations  in  the 
Assumed  Emittance  Values  of  the  Source  and/or  Substrate  Materials 

1.  Deviations  of  source  and  deposit  emittances  taken  together  - 

Temperature  distribution  for  the  deposition  of  yttrium  and  yttria  are 

calculated  for  various  emittance  values  of  the  condensate  and  the  evaporant  and 
are  tabulated  in  tables  VII  &  VIII.  As  the  deviation  in  emittances  increases 
the  errors  in  calculated  temperatures  increase  at  source  temperature  2200°K 
for  the  deposition  of  yttrium.  But  for  a  source  temperature  of  2400°K, 
and  hence  a  much  larger  evaporation  and  deposition  rates,  there  is  negligible 
errorless  than  1%,  even  for  the  largest  possible  error  in  emittances 
considered.  For  the  deposition  of  yttria,  the  errors  are  in  a  much  closer  range 
of  2  to  5%. 

2.  Deviations  of  substrate  emittance  values  only  -  Temperatures  on 
the  substrate  are  calculated  for  buffed  stainless  steel  substrate  and  rolled 


stainless  steel  substrate  and  are  tabulated  in  table  IX.  The  rolled  substrate 
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has  1.5  to  2  times  higher  emittance  than  buffed  substrate  and  hence  gives 
rise  to  5  to  8%  lower  temperatures  than  the  latter.  The  temperature  of  the 
source  and  the  distance  between  the  source  and  the  substrate  do  not  seem  to 
influence  the  difference  considerably. 

3.  Deviations  in  Emittance  of  the  Surroundings  of  the  Substrate 

Emittance  values  of  the  surroundings  are  varied  from  0  to  1  and  the 
substrate  temperatures  calculated  are  given  in  table  X.  It  is  observed 
that  when  the  substrate  temperature  is  greater  than  450°K,  the  errors  in 
temperatures  due  to  errors  in  the  emittance  of  the  surroundings  are  very 
small.  As  the  temperatures  of  the  substrate  increase,  the  errors  become 
less.  Around  500°K,  the  error  is  only  2%  even  when  the  emittance  of  the 
surroundings  varies  from  0  to  1.  At  600°K,  it  is  1%  and  at  700°K,  it  is 
<  1%. 

D.  Relative  Contributions  of  Various  Energy  Sources  in  Heating  of  the  Substrate 

The  various  energy  sources  contributing  to  the  heating  of  the  substrate 
are  the  radiant  energy  from  the  source,  the  energy  due  to  condensation  of 
the  vapors,  and  that  due  to  the  heat  of  reaction  (applies  to  yttria  only). 

Their  relative  contributions  are  calculated  for  various  experimental  condi¬ 
tions  for  the  deposition  of  yttrium  and  yttria  and  are  given  in  tables  XI  and 
XII.  They  are  a  function  of  the  material  being  deposited  and  the  temperature 
of  evaporation.  For  the  deposition  of  yttrium,  radiation  from  the  source 
contributes  twice  that  for  condensation  at  2100°K;  contributions  are  almost  equal 
at  2200°K.  At  2300°K,  the  contribution  of  condensation  is  twice  that  due 
to  radiation  and  at  2400°K,  the  former  contributes  about  3.5  times  the  latter. 
Thus,  we  see  that  there  is  a  reversal  of  relative  contributions  above  a  source 
temperature  of  2200°K.  The  relative  contribution  is  insensitive  to  the 


distance  between  the  source  and  the  substrate. 


For  the  deposition  of  yttria,  at  2000°K>,  radiation  contributes  10 
times  that  due  to  condensation  and  5  times  that  due  to  reaction.  As  the 
temperature  of  the  source  increased,  the  relative  contributions  vary.  Thus, 
at  2200°K,  radiation  contributes  two  times  the  condensation  and  same 
as  that  due  to  reaction.  Relative  contributions  of  radiation  and  condensation 
are  equal  and  reaction  contributes  two  times  either  one  at  2300°K.  Finally, 
at  2400°K,  condensation  contributes  two  times  and  the  reaction  contributes 
four  times  that  due  to  radiation. 

V.  SUMMARY  AND  CONCLUSIONS 

Temperature  and  thickness  variation  on  the  substrate  for  the  deposition 
of  yttrium  and  yttria  are  computed  for  various  experimental  conditions. 
Temperature  of  the  source,  diameter  of  the  source,  source  to  substrate  distance, 
material  of  the  substrate  and  the  thickness  of  the  substrate  are  the  important 
variables  in  these  experiments. 

There  are  three  sources  cf  heat  input  to  the  substrate,  viz.  radiation 
from  the  source,  condensation  and  reaction.  The  relative  contribution  of 
these  sources  is  calcula.ed  for  various  experimental  conditions.  Probable 
errors  in  temperatures  of  the  substrate  due  to  deviations  in  the  assumed 
emittance  values  of  the  substrate  material,  condensate  and  the  source  are 
estimated. 

One  of  the  stated  objectives  at  the  start  of  this  paper  was  to  use 
the  computations  to  focus  on  the  range  of  variables  to  be  explored  experimen¬ 
tally.  Let  us  now  set  some  criteria  for  the  production  of  a  desirable 
deposit  at  high  deposition  rates,  .0001-. 001"  per  idnute.  Such  a  deposit 
should  be  of  full  density  and  have  a  reasonably  uniform  microstructure 
(i.e.  variation  in  deposition  temperature  up  to  5%)  and  thickness  uniformity 
($  5%  from  center  to  edge). 
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From  the  calculations  and  discussion  above,  one  would  focus  on  the 
following  range  of  variables  to  be  explored  experimentally. 

1.  Temperature  of  the  source  -  2300  to  2600°K 

2.  Source  to  substrate  distance  -  25  to  35  cm. 

3.  Diameter  of  source  -  1.0"  to  1.5". 

4.  Substrate  material  -  buffed  stainless  steel  and  copper. 

5.  Substrate  thickness  -  .001  to  .005". 

6.  Substrate  size  (diameter)  -  8  to  20  cm. 

These  criteria  represent  compromises  between  the  possible  extremes. 

For  example,  one  can  get  great  thickness  uniformity  by  making  the  source-to- 
substrate  large  but  at  the  expense  of  deposition  rate  and  substrate  tempera¬ 
tures  which  may  be  too  low  to  produce  a  fully  dense  structure.  Furthermore, 
the  data  shows  that  errors  in  substrate  temperature  estimates  are  only 

-v 

slightly  influenced  by  variations  in  source  or  substrate  material  emittances 
by  operating  at  high  source  temperatures,  which  also  gives  high  deposition 
rates.  The  computations  also  illustrate  what  substrate  temperatures 
and  thickness  distributions  to  expect  in  case  certain  experimental  parameters 
are  fixed  or  beyond  one's  control.  In  such  cases,  one  might  resort  to  other 
methods  to  enhance  substrate  temperature  and  thickness  uniformity^  \  e.g. 
preheating  of  the  substrate^ beam-sweep  on  a  variable  residence  time/intensity 
pattern  in  a  rectangular  electron  beam  source  etc. 
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Appendix  -  1 


Vapor  Pressure  and  Evaporation  Rate  of  Yttrium 
Vapor  Pressure  Constants1*: 

log1Q  P  “  AT'1  +  B  logl0T  +  CT  +  DT2  +  E 

A  *  -25.2380  x  103 

B  “  -3.22193 

C  -  -1.18490  x  10"3 
D  =  1.64376  x  10"7 

E  *  23.4525 


• - 

Temperature 

°K 

Vapor  Pressure 
(mm) 

Evaporation  Rate 
g/cm2-sec. 

2,000 

0.305687  x  10'1 

0.375939  x  lO-3 

2,100 

0.926585  x  10_1 

0.111207  x  10'2 

2,200 

0.251396  x  10° 

0.294783  x  10-2 

2,300 

0.620095  x  10° 

0.711130  x  10'2 

2,400 

0.140913  x  101 

0.158198  x  10_1 

2.20 


Appendix  -  2 


Time  to  deposit  50  atomic  or  molecular  layers  on  the  substrate 

q 

Atom  Size  of  Yttrium  =  3.60  A° 

Equivalent  thickness  of  50  atoms  *  1.8  x  10  ^  cm. 


Time  to  achieve  50  atomic  layers 

1.8  x  10~6 _ 

condensation  rate 


sec. 


For  example 


At  source  temperature  of  22006K, 


H  ■  10  cm,  Delta 


Time 


1.8  x  10 


-6 


1.07  x  10 


-5 


0.168  sec. 


Temperature  of  the  Source  °K 


2200 

2300 

Thickness 

Thickness 

H(cm) 

Delta  (cm) 

rate 

cm/sec. 

Time 

(sec) 

rate 

cm/sec. 

Time 

(sec) 

10 

0 

1.07  x  10"5 

0.17 

2.58  x  10-5 

0.07 

15 

0 

1.15  x  10-5 

0.16 

0 

0.171  ::  10-5 

1.05 

0.413  x  10-5 

0.44 

25 

6 

0.153  10-5 

1.18 

0.369  x  10-5 

0.49 

50 

0 

0.427  x  10-6 

4.2 

0.103  x  10-5 

1.75 
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Heat  of  Condensation  at  298°K  for  Yttrium 


Heat  of  Sublimation  at  298°K  =  101.5  +  2.0  K  cal /mole. 

F  =  (Hm  -  H„.„)  +  heat  of  sublimation  at  298°K. 

T  298  vapor 


Temperature 

°K 

/U  _  U  N  ^-2 

'nT  298'vapor 

cal/mole 

F 

cal/mole 

1700 

7559 

109059 

1900 

8590 

110090 

2000 

9111 

110611 

2100 

9641* 

111141 

2200 

10172 

111672 

2300 

10720* 

112220 

2400 

11267 

112767 

2600 

12410 

113910 

Interpolated 
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TABLE  XIII 

Data  Used  In  These  Computations 


Symbol 


Parameter 


Value 


Reference 


M 

Molecular  weight  of  yttrium 

88.905 

13 

P 

Density  of  yttrium 

4.45  g/cm3 

13 

a 

Stefan-Boltzmann  constant 

1.356  x  10“12  cal/°K4-sec 

14 

h  or  ed 

Total  hemispherical  emmit- 
tance  of  Y  assumed  to  be 
the  same  as  spectral 
emittance  of  yttrium 
(X-0.65P) 

0.35 

15 

ed 

Total  hemispherical  emit- 
cance  of  Y2O3  assumed  to  be 
the  same  as  spectral  emit¬ 
tance  of  oxidized  yttrium 

0.60 

15 

e 

s 

Total  hemispherical  emic- 
tance  of  copper 

{0.01+(T-273)x9.09xl0-5} 

15 

e 

S 

Total  hemispherical  emit¬ 
tance  of  buffed  stainless 
steel 

(0.066+0. 267xTxl0_3) 

16 

e 

s 

Total  hemispherical  emit¬ 
tance  of  rolled  stainless 

(0.4025  +  10_4  x  T) 

16 

steel 


Specific  heat  of  copper 

(0.085  +  0.236  x  10_4xT) 

17 

Specific  heat  of  stain¬ 
less  steel 

(0. 120+0. 300x10_4xT)  cal/g. 

16 

Specific  heat  of  yttrium 

(5 . 59+1 . 9x10'3xT+0 . 29xl05/T2) 

17 

cal/deg-mole. 

Specific  heat  of  yttria 

(29.60+1.20xl0'3xT-4.78xl05/T2) 

17 

cal/deg-mole 

Heat  of  reaction  of 

2  Y(s)+  3/2  °2(g)  ■"  Y2°3(s) 

-AH“  =455.4  +  1.5 

K  cal/mole. 

18 

Density  of  copper 

8.96  g/cm3 

13 

Denisty  of  stainless  steel 

7.90  g/cm3 

19 
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2.  Synthesis  of  compounds  YjO^,  TiC,  Si^N^  by  reactive  evaporation 
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containing,  Gd,  La  and  Y  and  alloys  based  on  composition  Ni^Al!.^" 
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